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METHOD OF ANALYZING ELECTROMAGNETIC INTERFERENCE 



Background of the Invention 

1. Field of the Invention 

5 The present invention relates to a method of analyzing 

electromagnetic interference (EMI) (hereinafter often 
referred to as an "EMI analysis method") , and more particularly , 
to a method of analyzing electromagnetic interference arising 
in a large-scale, high-speed LSI (large-scale integrated 

10 circuit) by means of high-speed, highly-accurate logic 
simulation . 

2 . Description of the related Art 

LSIs find a broadening range of applications, from 
communications devices, such as cellular phones, to general 
15 household products, toys, and automobiles, as well as 

applications in the field of computers. Electromagnetic 
interference arising in such a product induces radio 
interference noise in a receiver, such as a TV set or a radio, 
or faulty operations of another system. In order to prevent 
20 this problem, the entirety of a product is shielded, or filters 
are provided in a product. With a view towards preventing an 
increase in the number of components and cost and difficulty 
in preventing occurrence of electromagnetic interference in 
a product, strong demand exists for suppression of noise in 
25 an LSI itself. 
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Under such a situation, an LSI is ranked as a key device 
for any product which contains an LSI. Demand exists for an 
larger-scale, high-speed LSI for ensuring competitiveness of 
a product. In a situation in which the cycle of product 
5 development becomes shorter, design-automation of an LSI is 
indispensable for satisfying the demand. There is growing 
necessity for adopting synchronous circuit as a condition for 
introducing a state-of-the-art design-automation method. In 
a case where all circuits of a large-scale, high-speed LSI 
10 operate synchronously with a reference clock signal, 

instantaneously-changing current becomes very large and 
induces an increase in electromagnetic interference. 

The present invention relates to a simulation method 
which enables evaluation of EMI indispensable for reducing 
15 electromagnetic interference while maintaining a tendency 
toward a larger-scale, higher-speed LSI. 

Noise imposed on another device by an LSI is roughly 
classified into two types; radiation noise, and conduction 
noise. Radiation noise emanated directly from an LSI includes 
20 noise emitted from internal wires of an LSI . However, internal 
wires do not act as an antenna of large size. As a matter of 
course, it is considered that the noise emitted directly from 
an LSI will pose a problem in the future, in association with 
an improvement in the operation frequency of an LSI . However , 
25 as of now, the noise emitted from the inside of an LSI is 
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considered trivial . 

In contrast, conduction noise affects another device 
mounted on a printed board, by way of direct interconnections, 
5 such as wires of an LSI or routings provided on a printed wiring 
board. Noise is emitted from such interconnections while the 
interconnections are act as the source of origination or as 
an antenna. The antenna constituted of the interconnections 
is much larger than that constituted by internal wires of an 
10 LSI and is a dominant element in terms of electromagnetic 
emission . 

A power line and a signal line can act as paths along 
which conduction noise developing in an LSI travels. In 
consideration of an electromagnetic field in the vicinity of 

15 an LSI, noise which results from variation in an electric 
current of a power source being emitted from a power line 
serving as an antenna is considered to be dominant. There may 
be a case where ringing and overshoot phenomena stemming from 
variation in a signal pose problems. However, there more 

20 frequently arises a case where variation in an internal power 
level of an LSI propagates as a signal waveform, to thereby 
present a problem. Noise emitted from a power line or a signal 
line is considered to have a strong correlation with variation 
in the electric current of a power source (hereinafter referred 

25 to as a "source current") . 
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A source current of a CMOS circuit will now be described 
by reference to a simple inverter circuit. In a case where 
variation arises in a voltage applied to an inverter circuit, 
there flows a load capacity charge/discharge current, which 

5 is the primary source current of the CMOS circuit . In addition , 
a short circuit current flows together with the load capacity 
char ge/dis charge current. In design of such a CMOS circuit, 
all circuits of an LSI are synchronized in accordance with 
constraints on the use of a design-automation tool . As a result 

10 of all circuits being synchronized, all circuits of the LSI 
operate simultaneously, and a peak current arises in a power 
source in synchronism with a reference clock signal . Further, 
in order to increase operating speed, or shorten a cycle, of 
the LSI, the capacity of a transistor is increased so as to 

15 enable a charging/discharging operation to.be completed within 
a short period of time. Eventually, a peak current increases. 
As a matter of course, the total source current of an LSI is 
increased when the level of an LSI is increased. Thus, the 
peak current of the power source is increased, thereby inducing 

20 occurrence of an abrupt change in a source current. Such an 
abrupt change induces an increase in higher harmonic components , 
thereby resulting in an increase in electromagnetic 
interference . 

Highly-precise simulation of change in a source current, 
25 which may be said to primarily account for electromagnetic 
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interference, is considered to be effective in evaluation of 
electromagnetic interference arising in an LSI. 

A current simulation method for effecting 
transistor- level current analysis, as will be described below, 
5 has conventionally been employed. 

FIG. 15 is a block diagram showing the flow of processing 
operations pertaining to a conventional transistor-level EMI 
analysis method. According to this method, on the basis of 
layout information pertaining to an LSI which is to be analyzed 
10 through use of a transistor -level current analysis method, 
there is performed layout parameter extraction (hereinafter 
referred to simply as an "LPE") processing 03. Subsequently, 
there is performed circuit simulation 06 regarding a 
switch-level netlist; source-of -current modeling 08; a power 
15 line LPE step O10; transient analysis simulation 012; and FFT 
processing 014 . 

Processing pertaining to each of the foregoing 
processing steps will now be described by reference to FIG. 
15 . 

20 First, in step 03 are input layout data Ol pertaining 

to a semiconductor integrated circuit to be subjected to EMI 
analysis; parameters of elements, such as transistor elements 
or various parasitic wiring elements (e.g., resistors and 
capacitors) ; and an LPE rule 02 for defining a form in which 

25 extracted layout parameters are to be output. In accordance 
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with the LPE rule 02 , parameters of the respective elements 
included in the layout data Ol are calculated, whereby a netlist 
04 is produced. In step 03 , parasitic elements of a power 
source (and the ground) are not objects of extraction. 
5 In step 06 are input the netlist 04 prepared in step 03 

and a test pattern 05 for causing a circuit, which serves as 
an object of analysis, to replicate a desired logic operation. 
There are calculated a load capacity charge/discharge current 
and a short circuit current, which correspond to the operating 
10 state of an internal circuit, thereby producing current 

waveform information 07 concerning the waveform of an electric 
current of a transistor . The first operation of the processing 
pertaining to step 06 is effected on the assumption that the 
potential of a power source (and that of ground) is a 
15 variation-free, ideal potential. 

In step 08 is entered the current waveform information 
07 concerning a transistor prepared in step 06. The thus- 
entered current waveform information 07 is modeled into a mode 
which can be applied to subsequent step 012, wherewith current 
20 source element model information 09 is prepared. In order to 
alleviate a load which would be imposed on subsequent step Ol2 , 
a function circuit block consisting of a plurality of 
transistors is usually modeled as a single current- source 
element . 

25 Processing pertaining to step O10 differs from 
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processing pertaining to step 03, only in that rather than 
parameters of transistor elements and those of various 
parasitic wiring elements, parameters of parasitic elements 
of a power source and those of a ground wire (e.g. , resistors, 
5 decoupling capacitance, and like elements) are taken as objects 
of extraction. Hence, repeated explanation is omitted. In 
step O10, a power source (and ground) wiring netlist Oil is 
produced . 

In step 012 are entered the current source element model 

10 information 09 prepared in step 08 , the power source (and 
ground) wiring netlist Oil prepared in step 010, and impedance 
016 of a wire or a lead frame (including, resistance, 
capacitance, and inductance) . Through analysis of these input 
data carried out by a transient analysis simulator typified 

15 by SPICE, fluctuations in line voltage of a circuit to be 
analyzed are calculated. Thus, there is produced a line 
voltage drop result 017 concerning the thus-calculated 
fluctuations in line voltage. 

Subsequently, processing pertaining to step 06 is 

20 performed again. In contrast with the first operation of the 
processing pertaining to step 06 having been effected on the 
assumption that the potential of the power source (and the 
ground) is a fluctuation-free, ideal potential, the line 
voltage drop result 017 prepared in step 012 is entered. The 

25 current waveform information 07 concerning a transistor is 
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prepared again in consideration of fluctuations in line voltage. 
Similarly, processing pertaining to steps 08 and 012 is 
repeated . 

Processing pertaining to steps 06, 08 , and 012 is 
5 effected several times in a looped manner, wherewith there is 
produced a current waveform result 013 which highly-accurately 
duplicates fluctuations in line voltage. 

In step 014 , the current waveform result 013 prepared 
in step 012 is entered and subjected to fast Fourier 
10 transformation (hereinafter abbreviated FFT) , to thereby 

enable frequency spectrum analysis. There is obtained an EMI 
analysis result Ol5. 

In the conventional example, the precision of 
verification varies greatly according to combination of the 
15 LPE processing 03, the power line LPE processing 010, and the 
source current modeling processing 08 . However, a certain 
level of accuracy of analysis can be expected. A transient 
analysis simulator typified by SPICE is used for 
trans is tor- level analysis of an electric current. Hence, a 
20 limitation is imposed on the level of a circuit to be analyzed, 
and an enormous amount of processing time is required. The 
level of a semiconductor integrated circuit has increased 
recently, and establishment of an EMI analysis method which 
enables high-speed analysis of an electric current on a level 
25 larger than a transistor level is desired. 
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A gate-level current analysis method has conventionally 
been proposed as a current analysis method which can be made 
faster. This gate-level current analysis method is used for 
analyzing power consumption . One example of a gate-level 
5 current analysis method is EMI -noise analysis which is to be 
effected in an ASIC design environment. The method is 
described in "EMI-Noise Analysis Under ASIC Design 
Environment" (ISPD&99, pp. 16 through 21) . According to this 
technique, an event is acquired from the result of a gate- 

10 level simulation using a test vector, and the waveform of an 
electric current is estimated. The frequency of the thus- 
estimated current waveform is analyzed through fast Fourier 
transformation (FFT) . More specifically, as shown in FIG. 16, 
a logic simulation 104 is effected on the basis of a netlist 

15 101 and a test vector 102, wherewith event information 105 is 
calculated. On the basis of the thus -calculated event 
information 105 and waveform information 103 obtained at the 
time of toggling, processing pertaining to a current waveform 
calculation section 107 is executed, to thereby produce a 

20 current waveform calculation result 108. This current 

waveform calculation result 108 is subjected to FFT processing 
109, to thereby produce a frequency characteristic 110. The 
EMI-noise analysis method can effect an EMI analysis operation 
faster than that performed according to the conventional gate 

25 - level EMI analysis method. However, use of a test vector 



9 



still involves consumption of much execution time. Therefore, 
the processing speed achieved by the EMI -noise analysis method 
is not sufficiently high, and demand still exists for an 
increase in processing speed of the EMI-noise analysis method. 
The EMI-noise analysis method also encounters a problem of an 
analysis result being dependent on the pattern of an employed 
test vector. 

As mentioned above, the conventional example using the 
transistor-level current analysis method can be expected to 
yield a certain level of accuracy. However, a transient 
analysis simulator typified by SPICE is used for 
transistor-level current analysis. A limitation is imposed 
on the level of a circuit to be analyzed, and an enormous amount 
of processing time is required. The level of a semiconductor 
integrated circuit has increased recently, and there is desired 
establishment of an EMI analysis method which enables 
high-speed analysis of an electric current at a scale larger 
than the scale which can be analyzed by a transistor-level 
simulator . 

The gate-level simulation using a test vector has also 
been proposed. However, the example conventional gate-level 
simulation technique encounters difficulty in increasing the 
speed of analysis . Since the gate-level simulation technique 
employs a test vector, an analysis result is dependent on the 
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employed test pattern. 

Summary of the Invention 
The present invention has been conceived to solve the 
drawbacks of the conventional methods and is aimed at 
evaluating electromagnetic interference developing in an LSI 
through a simulation by means of high-speed, highly-accurate 
analysis of a power-supply current. 

To this end, the present invention provides a method of 
analyzing electromagnetic interference (an EMI analysis 
method) . In contrast with a known dynamic gate-level 
simulation method, the EMI analysis method enables estimation 
of EMI noise, by means of calculating signal propagation of 
each node through use of the signal propagation probability 
technique, and calculating variation time of each node through 
use of "the Static timing analysis technique". In short, the 
present invention is characterized in calculating a frequency 
characteristic from the relationship between toggle 
probability of each node and delay in each node. 

More specifically, the present invention provides a 
method of analyzing electromagnetic interference generating 
in an LSI, comprising: 

a correction step of correcting the amplitude of a 
current estimation waveform in each node which has been 
previously prepared for each change in each node, in accordance 
with the probability of variation in each node; 
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an addition step of adding current waveforms of all nodes 
together within a period of time corresponding to one cycle, 
provided that the thus-corrected current waveform appears at 
a time a signal arrives at each node; and 
5 a frequency analysis step of analyzing the frequency of 

the current waveform calculated in the addition step. 

The probability of variation in each node is calculated 
through use of the signal propagation probability technique. 
Further, the time at which a signal arrives at each node is 
10 calculated through use of the static timing analysis technique . 
The time at which a signal arrives is defined so as to fall 
within the range between the maximum time and the minimum time, 
in accordance with an average time, the maximum time, the 
minimum time, or a predetermined distribution such as a normal 
15 distribution. 

Preferably, in the correction step, the amplitude of a 
current estimation waveform, which has been prepared for each 
change in each node, is corrected in accordance with the 
probability of variation in each node and a distribution with 
20 respect to time (hereinafter called "chronological 
distribution") . 

More specifically, the probability of change in each node 
is calculated through use of the signal propagation probability 
technique, and the chronological distribution at which a signal 
25 arrives at each node is calculated through use of the static 
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timing analysis technique. 

Preferably, each node has a plurality of paths, and a 
current waveform is calculated in consideration of a case where 
each of the paths has a unique probability of change and signal 
5 arrival time . 

Further, the present invention provides a method of 
analyzing electromagnetic interference developing in an LSI, 
the method comprising: 

a waveform formation step of forming a current estimation 
10 waveform which has been prepared for each change in each node, 
as if the waveform randomly arises within a plurality of 
predetermined cycles, in accordance with the probability of 
change in each node and a time at which a signal arrives at 
each node; 

15 adding the thus -prepared current estimation waveforms 

of all nodes, to thereby derive a current waveform; and 

analyzing the frequency of the current waveform, thereby 
determining a noise characteristic of EMI . 

In other words , under the EMI analysis method according 

20 to the present invention, the probability of change in each 
node is calculated through use of the signal propagation 
probability technique, and the result of calculation is stored 
as a probability at which a signal randomly changes . Further, 
a time at which a signal arrives at each node is calculated 

25 through use of the static timing analysis technique. 
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Moreover, the present invention provides a method of 
analyzing electromagnetic interference developing in an LSI, 
the method comprising: 

a waveform formation step of forming a current estimation 
5 waveform which has been prepared for each change in each node, 
as if the waveform randomly arises within a plurality of 
predetermined cycles, in accordance with the probability of 
change in each node and chronological distribution 
probability; 

10 adding the thus-prepared current estimation waveforms 

of all nodes, to thereby derive a current waveform; . and 

analyzing the frequency of the current waveform, thereby 
determining a noise characteristic of EMI. 

More specifically, the probability of change in each node 

15 is calculated through use of the signal propagation probability 
technique, and the result of calculation is stored as the 
probability of a signal changing randomly. A chronological 
distribution at which a signal arrives at each node is 
calculated through use of the static timing analysis technique. 

20 Preferably, each node has a plurality of paths, and a 

current waveform is calculated in consideration of a case where 
each of the paths has a unique probability of change and signal 
arrival time . 

Brief Description of the Drawings 
25 FIG. 1 is a schematic illustration for describing the 
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concept of the present invention; 

FIG. 2 is a block diagram showing a portion of a circuit 
used in a first embodiment of the present invention; 

FIGS. 3A and 3B are graphs showing waveforms of signals 
5 arriving at the respective nodes of the cell shown in FIG. 2; 

FIG. 4 is a block diagram showing the processing of a 
frequency characteristic calculation block according to a 
first embodiment of the present invention; 

FIGS . 5A to 5D are illustrations for describing a 
10 processing image according to the first embodiment of the 
present invention; 

FIG. 6 is a flowchart of processing of a current waveform 
calculation according to the first embodiment; 

FIG. 7 is a block diagram showing the processing of a 
15 frequency characteristic calculation block according to a 
second embodiment of the present invention; 

FIGS. 8A through 8D are illustrations for describing a 
processing image pertaining to the second embodiment; 

FIG. 9 is a flowchart of processing of a current waveform 
20 calculation according to the second embodiment; 

FIG. 10 is a block diagram showing a portion of a circuit 
used in a third embodiment of the present invention; 

FIG. 11 is a graph showing signal waveforms of each path; 
FIG. 12 is a block diagram showing a frequency 
25 characteristic calculation block according to a third 
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embodiment of the present invention; 

FIGS. 13A through 13C are illustrations showing a 
processing pertaining image to the third embodiment; 

FIG. 14 is a flowchart of current waveform calculation 
5 processing according to the third embodiment; 

FIG. 15 is a flowchart for describing a known EMI analysis 
method; and 

FIG. 16 is a flowchart for describing a method of 
analyzing EMI dynamic at gate-level . 
10 Description of the preferred Embodiments 

An electromagnetic interference analysis method 
according to preferred embodiments of the present invention 
will now be described by reference to the accompanying drawings . 
As shown in FIG. 1, an EMI analysis method according to the 
15 present invention is characterized in: 

calculating the transition probability of a node from 
a netlist 1 and a transition probability 2, through use of a 
propagation probability method, and calculating a static delay 
4 through use of a static delay analysis method, to thereby 
20 derive a calculated probability/delay 5 of a node; 

estimating the waveform 6 of an electric current on the 
basis of the probability/delay 5 and information 3 concerning 
the waveform of an electric signal at the time of toggling, 
to thereby derive a current waveform estimation result 7; and 
25 subjecting the current waveform estimation result 7 to 



16 



a fast Fourier transformation 8 (hereinafter called an "FFT") , 
thereby determining a frequency characteristic 9 of the 
waveform. 
(First Embodiment) 

A method of analyzing electromagnetic interference 
according to a first embodiment of the present invention will 
be described hereinbelow. As can be seen from a schematic 
diagram shown in FIG. 1, under an EMI analysis method according 
to the present embodiment, 

the quantity of electromagnetic interference developing in an 
LSI is to be analyzed on the basis of a transient probability 
and static delay propagation data, provided that a waveform 
shown in FIG. 3A appears in a node A of a flip-flop (FF) cell 
and a waveform shown in FIG. 3B appears in a node B of the FF 
cell (where FIG. 3B is an enlarged view of about 1.5 cycles 
of the signal designated by braces in FIG. 3A) when a clock 
signal CLK is input to a circuit shown in FIG. 2. Here, the 
transition probability of a node is calculated from a 
previously-prepared netlist 1 and a transition probability 2. 
Further, a static delay 4 in a current estimation waveform per 
change is calculated. The amplitude of a current waveform is 
corrected in consideration of information 3 concerning the 
waveform of an electric current arising at the time of a 
predetermined toggling operation. Provided that the 
corrected current waveform arises at a time at which a signal 
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arrives at the respective node, the current waveforms which 
appear at all nodes during a period of time corresponding to 
one cycle are added (the current waveform estimation processing 
6) . The current waveform estimation result 7 determined 
through addition is subjected to the FFT processing 8, thereby 
determining the frequency characteristic 9 of EMI components 
of a circuit to be analyzed. 

FIG. 4 is a block diagram for describing the overall flow 
of processing of the EMI analysis method according to the 
present embodiment. FIGS. 5A through 5D are illustrations 
showing the principle underlying the processing. In a netlist 
401 , a circuit which is an object of EMI analysis is represented 
as circuit data. Delay information 405 concerning each node 
is formed from the netlist 401 through static delay calculation 

403 (see FIG. 5A) . Transition probability information 406 
concerning each node is formed from the netlist 401 and input 
transition probability 402, through propagation probability 

404 (see FIG. 5B) . In consideration of a triangular waveform 
whose area corresponds to the quantity of electric current 
derived by means of multiplying current waveform information 
by probability information, an average current waveform 4 09 
is formed by average current waveform calculation means 408 
from element current waveform information 407 concerning each 
node (see FIG. 5C) and the delay information 405. The 
thus -determined average current waveform 4 09 is taken as 
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average current waveform information (see FIG. 5D) . The 
average current waveform information is subjected to FFT 
processing 410, thereby deriving frequency characteristic 
information 411. 

FIG. 6 shows a flowchart of processing of the average 
current waveform calculation means 408. The average current 
waveform calculation means 408 reads element current waveform 
information from a table (step 1250) and performs a current 
waveform calculation loop (step 1251) . The base of a 
triangular waveform of an instance to be processed is extracted 
from an output slew (step 1252) . The area of the triangular 

h 

waveform is taken as being derived by means of multiplying W x — 

by transition probability per cycle, and I is taken as the value 
of the area of the triangular waveform. The height "h" of the 
triangular waveform is calculated from transition probability 

per 2x1 cycle (step 1253) , wherein "I" denotes the quantity 
^ Wxl 

of electric current flowing in a cell of an event which is an 
object of processing. This processing corresponds to 
processing performed by a triangular waveform shaping section. 
Until variable «x" changes from 0 to W/2 , h(c, i) expressed 
(by Eq.5) is added to I (t+x) and I (t-x) . Further, At is added 
to variable "x" (steps 1254 and 1255). Here, I (t+x) denotes 
total electric current flowing through all the cells at time 
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t+x, and I (t-x) denotes total electric current flowing through 
all the cells at time t-x. 

The frequency characteristic of a circuit to be analyzed 
can be determined in the manner as mentioned previously, and 
5 a designer can analyze EMI which would arise in a circuit of 
interest . 

According to the EMI analysis method, a current waveform 
is modeled through an averaging operation, on the basis of 
static delay information and propagation probability 
10 information. The thus-obtained model is subjected to FFT 
processing, thereby analyzing EMI of a circuit.. The EMI 
analysis method can analyze EMI components within a shorter 
period of time than can a known gate-level dynamic analysis 
method. 

15 i n a case where performance of EMI analysis for each path 

of a circuit is desired, static delay information concerning 
each path is given. 

In the present embodiment, electric currents of all nodes 
in a circuit to be analyzed are added. However, so long as 

20 the number of nodes whose electric currents are to be added 
is controlled in accordance with the magnitude of an electric 
current or the frequency of probability, processing time can 
be shortened further. 
(Second Embodiment) 

25 Next will be described an EMI analysis method according 
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to a second embodiment of the present invention . As 
represented by a flowchart shown in FIG. 7, the present EMI 
analysis method is characterized in employing random current 
waveform estimation means 708 in lieu of the average current 
5 waveform calculation means 4 08 employed in the first embodiment, 
and utilizing random current waveform information in lieu of 
the average current waveform information. In other respects, 
the EMI analysis method according to the present embodiment 
is identical in configuration with that described in connection 

10 with the first embodiment. 

FIG. 7 a block diagram for describing the overall flow 
of processing of the EMI analysis method according to the 
present embodiment. FIGS. 8A through 8D are illustrations 
showing the principle underlying the processing. In a netlist 

15 701, a. circuit which is an object of EMI analysis is represented 
as circuit data. Delay information 705 concerning each node 
is formed from the netlist 701 through static delay calculation 

703 (see FIG. 8A) . Transition probability information 706 
concerning each node is formed from the netlist 701 and input 

20 transition probability 702, through propagation probability 

704 (see FIG. 8B) . On the basis of element current waveform 
information 707 concerning each node (FIG. 8C) and operating 
frequency information 712, random waveform estimation means 
708 produces random current waveform information 709 (see FIG. 

25 8D) within a plurality of predetermined cycles. The thus- 
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produced random current waveform information 709 is subjected 
to FFT processing 710, thereby deriving frequency 
characteristic information 711. 

FIG. 9 shows a flowchart of processing of the random 
5 current waveform estimation means 708 . The average current 
waveform estimation means 708 reads element current waveform 
information from a table (step 1280) and performs a current 
waveform calculation loop (step 1281) . The average current 
waveform estimation means 708 performs loop processing until 

10 valuable "y" changes from 1 to the value of a. frequency, (step 
1282) . The following processing is iterated until calculation 
of a current waveform is completed. A determination is made 
as to whether or not a random number is smaller than the value 
of probability (step 1283) . If a random number is smaller, 

15 the base of a triangular waveform of an instance to be processed 
is extracted from an output slew (step 1284) . At this time, 

the area of the triangular waveform is defined as W x-^- , and 

I is the value of the area of the triangular waveform. The 

height "h" of the triangular waveform is calculated by 2x ^r 

20 (step 1285) , wherein "I" denotes the quantity of electric 
current flowing in a cell of an event which is an object of 
processing. This processing corresponds to processing 
performed by a triangular waveform shaping section. Until 
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variable "x" changes f rom 0 to W/2 , h(c, i) expressed (byEq.5) 
is added to I (t+x) and I (t-x) . Further, At is added to variable 
"x" (steps 1286 and 1287) . Here, I (t+x) denotes total electric 
current flowing through all the cells at time t+x, and I (t-x) 
5 denotes total electric current flowing through all the cells 
at time t-x . 

The frequency characteristic of a circuit to be analyzed 
can be determined in the manner as mentioned previously, and 
a designer can analyze EMI which would arise in a circuit of 
10 interest. 

According to the EMI analysis method, a current waveform 
is modeled through a random current waveform operation, on the 
basis of static delay information and propagation probability 
information. The thus-obtained model is subjected to FFT 

15 processing, thereby analyzing EMI of a circuit. The EMI 

analysis method can analyze EMI components with high accuracy 
within a shorter period of time than a known gate-level dynamic 
analysis method. 

In a case where performance of EMI analysis for each path 

20 of a circuit is desired, static delay information concerning 
each path i s given . 

In the present embodiment, electric currents information 
of all nodes in a circuit to be analyzed are added. However, 
so long as the number of nodes whose electric currents are to 

25 be added is controlled in accordance with the magnitude of an 
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electric current or the frequency of probability, processing 
time can be shortened further. 
(Third Embodiment) 

An EMI analysis method according to a third embodiment 
5 of the present invention will now be described. In the 
previously-described first and second embodiments, delay 
information and probability information are prepared 
separately. Information is derived by means of multiplying 
waveform information which is obtained as element current 
10 waveform information, by probability information. The. 

thus -obtained information is added to a delay time of each node . 
In contrast, in the present embodiment, delay propagation 
probability information is formed from delay propagation 
probability information. Delay/ transition probability is 
15 calculated from the delay propagation probability information , 
and element waveform information is added to the thus- 
calculated del ay/ transit! on probability. 

In this way, more realistic current waveform information 
is calculated, and the result of current waveform calculation 
20 is subjected to FFT processing, thereby determining the 

frequency characteristic of an EMI component of a circuit to 
be analyzed. Thus, EMI of the circuit is analyzed. As can be 
seen from an enlarged view shown in FIG. 10, the present 
embodiment is directed particularly a case where a plurality 
25 of paths are provided in a composite cell . FIG. 11 shows delay 
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transition information concerning propagation of a signal in 
each of the paths of the composite cell shown in FIG. 10. FIG. 
11 shows delay transition probability information as one 
example. As can be seen from FIG. 11 , there is obtained node 
5 information concerning a plurality of paths, and mean current 
waveform information is formed from the node information. 

FIG. 12 is a block diagram for describing the overall 
flow of processing of the EMI analysis method according to the 
present embodiment. FIGS. 13A through 13C are illustrations 

10 showing the principle underlying the processing.. FIG. 14 is 
a flowchart of average current waveform calculation means used 
in the processing. In a netlist 901, a circuit which is an 
object of EMI analysis is represented as circuit data. 
Delay/transition probability 906 of each node is calculated 

15 from the netlist 901 and input transition probability 902, on 
the basis of delay/propagation probability 904 (see FIG. 13A) . 
Mean current waveform estimation means 908 produces mean 
current waveform information 909 (see FIG. 13C) , in 
consideration of a triangular waveform whose area is determined 

20 by the quantity of electric current, such that the 

del ay /transit! on probability 906 is multiplied by element 
current waveform information 907 (see FIG. 13B) . The 
thus-f ormed mean current waveform information 909 is subjected 
to FFT processing 910 within a time domain which is determined 

25 on the basis of operating frequency information 912, thereby 
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obtaining frequency characteristic information 911 . 

FIG. 14 shows a flowchart of processing of the average 
current waveform calculation means . The average current 

5 waveform calculation means reads element current waveform 
information from a table (step 1310) and performs a current 
waveform calculation loop (step 1311) . The following 
processing is iterated until calculation of a current waveform 
is completed. The delay/ transition probability 906 

10 calculated from delay information and transition probability 
information is multiplied by element current waveform 
information 907 (see FIG. 13B) (step 1312) . In consideration 
of a triangular waveform whose area is determined by the 
quantity of electric current, average electric current 

15 waveform estimation means 908 adds the result of multiplication 
as mean current, thereby deriving average current waveform 
information 909 . The average current waveform information 909 
is subjected to FFT processing 910, thereby determining 
frequency characteristic information 911. 

20 

The frequency characteristic of a circuit to be analyzed 
can be determined in the manner as mentioned previously, and 
a designer can analyze EMI, which would arise in a circuit of 
interest. 

25 
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According to the EMI analysis method, delay propagation 
probability information is formed from static delay 
information and propagation probability information, and 
average current waveform information is formed from the delay 
5 propagation probability information. The thus-obtained 
average current waveform information is subjected to FFT 
processing, thereby enabling highly-accurate EMI analysis. 
The EMI analysis method can analyze EMI components within a 
shorter period of time than can a known gate-level dynamic 
10 analysis method. 



In addition to a distribution taking into consideration 
a path, a temperature/process/voltage distribution may be 
conceived as the delay/ transition probability information, 
shown in FIG. 13A. 



In the foregoing embodiments, FFT processing has been 
used for analyzing a frequency. However, the present 
invention is not limited to FFT processing. Needless to say, 
20 another processing method, such as ordinary Fourier 
transformation, may alternatively be employed. 



The present invention can materialize evaluation of EMI 
25 developing in an LST through a simulation, by means of 
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highly-accurately analyzing through static processing, 
variation in power-supply current which may be said to 
primarily account for EMI. Further, in contrast with dynamic 
analysis of EMI which is embodied by a gate-level simulation 
or a like simulation, the present EMI analysis method can 
prevent an increase in processing time. 
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What is claimed is: 

1. A method of analyzing electromagnetic interference 
developing in an LSI, comprising: 

a correction step of correcting the amplitude of a 
current estimation waveform in each node which has been 
previously prepared for each change in each node, in accordance 
with the probability of variation in each node; 

an addition step of adding current waveforms of all nodes 
together within a period of time corresponding to one cycle, 
provided that the thus -corrected current waveform appears at 
a time a signal arrives at each node; and 

a frequency analysis step of analyzing the frequency of 
the current waveform calculated in the addition step. 

2 . The method of analyzing electromagnetic interference 
developing in an LSI according to claim 1, wherein the 
correction step includes a step of correcting the amplitude 
of a current estimation waveform, which has been prepared for 
each change in each node, in accordance with the probability 
of variation in each node and a distribution with respect to 
time . 

3. The method of analyzing electromagnetic interference 
developing in an LSI according to claim 1, wherein each node 
has a plurality of signal transmission paths (hereinafter 
referred to simply as "paths") , and each of the current 
waveforms is calculated in consideration of a case where each 
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of the paths has a unique probability of change and signal 
arrival time . 

4 . The method of analyzing electromagnetic interference 
developing in an LSI according to claim 2, wherein each node 
5 has a plurality of paths, and each of the current waveform 
is calculated in consideration of a case where each of the paths 
has a unique probability of change and signal arrival time. 

5. A method of analyzing electromagnetic interference 
developing in an LSI , the method comprising: 
10 a waveform formation step of forming a current estimation 

waveform which has been prepared for each change in each node, 
as if the waveform randomly arises within a plurality of 
predetermined cycles, in accordance with the probability of 
change in each node and a time at which a signal arrives at 
15 each node; 

adding the thus-prepared current estimation waveforms 
of all nodes, to thereby derive a current waveform; and 

analyzing the frequency of the current waveform, thereby 
determining a noise characteristic of EMI . 
20 6. The method of analyzing electromagnetic interference 

developing in an LSI according to claim 5, wherein each node 
has a plurality of paths, and a current waveform is calculated 
in consideration of a case where each of the paths has a unique 
probability of change and signal arrival time. 
25 7 . A method of analyzing electromagnetic interference 
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developing in an LSI, the method comprising: 

a wavef orm formation step of forming a current estimation 
waveform which has been prepared for each change in each node, 
as if the waveform randomly arises within a plurality of 
5 predetermined cycles, in accordance with the probability of 
change in each node and a distribution probability of time; 

adding the thus-prepared current estimation waveforms 
of all nodes, to thereby derive a current waveform; and 

analyzing the frequency of the current waveform, thereby 
10. determining a noise characteristic of EMI. 

8 . The method of analyzing electromagnetic interference 
developing in an LSI according to claim 7, wherein each node 
has a plurality of paths, and a current waveform is calculated 
in consideration of a case where each of the paths has a unique 
15 probability of change and signal arrival time. 
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Abstract 

In contrast with a known dynamic gate-level simulation method, 
a method of analyzing electromagnetic interference (an EMI 
analysis method) according to the present invention enables 
estimation of EMI noise, by means of calculating signal 
propagation of each node through use of the signal propagation 
probability technique, and calculating variation time of each 
node through use of "the Static timing analysis technique". 
In short, the present invention is characterized in calculating 
a frequency characteristic from the relationship between 
toggle probability of each node and delay in each node. 
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